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شبیه‌سازی و بهینه‌سازی حل عددی در تحلیل توزیع دما و تنشهای پسماند حرارتی در پوشش دولایه 
سد حرارتی 8:77:007/8۷57,]به کمک حل المان محدود* 
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نسیم نایب پاشایی حسین وفایی نژاد سید حسین سیدین محمدرضا اپوطالبی 


سید محمد مهدی هادوی " حسین سرپولکی(* 
توزیع دما و تنش‌های پسماند در پوشش سدحرارتی سه لایه ( 022۲)/8۲52/[۷)۲۸۵۲) تحت یک چرحه حرارتی ت رکیبی و واقعی به 
کمک نرم افزار تجاری اباکوس شبیه‌سازی شد. از تکنیک کاهش زمان حل به روش مفیأس‌دهی جرمی برای جلوگیری از اعوجاج بیش از حد 
مشش‌بندی, کاهش خحطای عددی و واکرایی استفاده گردید. نتایج نشان داد که استفاده از روش افزایش جرم المان‌ها و جابجای یگره‌ها بصورت 
تطبیق‌پذیر باعث افزایش سرعت حل می‌شود که دلی لآن بزرگتر شدن نموهای حلی در فرایند است. نتایج شبیه‌سازی نشان داد که تنش‌های 
پسماند در منطقه پوشش محافظ و نه زیرلایه تم رکز یافته است و این به نویه حود باعث بهبود و افزایش عمر این محصول محافظ می‌باشد. 


تحلیل عددی نشان داد که بیشترین میزان اعوجاج بصورت عمده در ناحیه پوث ۷ ای سرامیکی متم رکز شده است. 


واژه‌های کلیدی پوشش سد حرارتی» تنش پسماند. شبیه سازی المان محدود. شوک حرارتی 


ص مع‌ووما5 امعم له صمتاب‌طاتتعزنا هه مره 1 ۵ جمتامابنه )صعصعا۴ مانصز۳ 
۸ ۱,212 ما۵۵ ۵۲۲۱۵۲ امصصمط 1 م۱ 


نام نویسندگان به انگلیسی ارسال گردد 


و 


مد امصصعط؟ رصمتتمابنه تمعصوماه ماتصز۲ روعومم‌تاو لمبل‌توم؟ رعصتلعمی تمتنده تقحصعط 1 ۷۷۵۲۵8 167 


*نسخه‌ی نخست مقاله در تاریخ و نسخه ی پایانی آن در تاریخ به دفتر نشریه رسیده است. 


(۱) نویسنده مسئول: دانشجوی دکتری, دانشکده مهندسی مواد و متالورژی, دانشگاه علم و صنعت ایران 
(۲) دانشجوی دکتری. دانشکده مهندسی مواد و متالورژی, دانشگاه علم و صنعت ایران 

(۳) دانشیار دانشکده مهندسی مواد و متالورژی دانشگاه علم و صنعت ایران 

(۶) استاد. دانشکده مهندسی مواد و متالورژی, دانشگاه علم و صنعت ایران 

(۵) دانشیار دانشکده مهندسی مواد و متالورژی. دانشگاه مالک اشتر 


(1) استاد دانشکده مهندسی مواد و متالورژی, دانشگاه علم و صنعت ایران 
:1۳01 


مقدمه 
سوپرآلیاژهای مورد استفاده در اجزای توربین‌های دما 
بالاء به دلایلی همچون اکسیداسیون و گازهای خورنده 
که در شرایط کاری و در نتیجه احتراق آزاد می‌شوند 
در معرض شرایط کاری دشوار قرار دارند [1-1]. 
ترکینت شیمیانی این عویر الیازها به کونه‌ای اشست کنه 
می‌توانند استحکام خوب و مقاومت خزش در دمای بالا 
را تافین کل اما به تهاین حفاطت مخیظین کافی: ندار نت 
ارات ادا توش یاقب ام ای سوه 
زیرلایه از تخریب در شرایط کاری سخت ضروری 
است [۵01]. تلاش‌های متعددی در دهه‌های گذشته 
منجر به توسعه پوشش‌های حفاظتی برای اجزای توربین 
شده است [/۵0۷]. 

پوشش‌های سد حرارتی وظیفه کاهش دمای اجزای 
موجود در بخش‌های داغ توربین‌های گازی و در نتیجه 
افزایش طول عمر آن‌ها را بر عهده دارند [۸4]. کاربرد 
ان تفش هرت کتا مهن تی احفتال را ادن 
کاسته شدن از دمای زیرلایه می‌شود که بهبود مصرف 
سوخت. افزايش قدرت موتور و ماندگاری اجزای 
موتور را در پی خواهد داشت. تحقیقات نشان داده 
است که این پوشش‌ها می‌توانند مقاومت خوبی در 
برابر سایش و خوردگی ایجاد کنند [۱۷۱۰]. 

عن فک رشن هام بل کر آرزتی مدا تصتي 
تاثیر ترکیب شیمیایی لایه واسط نوع اکسید محافظ دما 
و تغییرات آن. روش اعمال پوشش عوامل محبطی 
مانند ترکیب شیمیایی و سرعت گاز نمک مذاب 
خورنده و فرسایش ناشی از گازهای داغْ است [۶-۱1]. 
وظیفه اصلی پوشش‌های سد حرارتی جلوگیری از 
انتقال جریان گرمایی به زیرلایه است. بنابراین» قابلیت 
عایق سازی حرارتی به عنوان یکی از مهمترین فاکتورها 
در ارزیابی عملکرد و توسعه کاربردی این پوشش‌ها به 
شمار می‌رود [۱۱۰۱۲]. 


پوشش‌های سد حرارتی زیرکونیای به طور جزئشی 
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پانشان ه با اسر( )اور فش ای ۱۳۳۵۲۲ 
(۷- 00 و لا[ ویا ترکیبی از این دو عنصر) به دلیل دارا 
بودن طیفی از خواص موردنیاز. بیشترین کاربرد را به 
عنوان پوشش‌های محافظ در دمای بالا دارند 
۸۱۳ در کنار مزایای فراوان این دسته از پوشش- 
های محافظ نقطه ضعف اصلی ۷۲۹2 دمای عملیاتی 
محدود (۱۲۰۰۹۵) برای کاربردهای طولانی مدت است 
[۱۶۰۱۵]. در دماهای بالات استحاله فازی فاز شبه پایدار 
تراگونان (0) به فاز تتراگونال و فاز مکعبی باعث مسی- 
شود که حین سرمایش فاز تتراگونال به فاز منوکلینیک 
تبدیل شود و میکروترک‌ها در پوشش تشکیل شوند. 
علاوه بر اين. با توجه به زینتر شدن پوشش در دماهای 
بالاه تخلخل پوشش کاهش می یابد در حالی که مدول 
یانگ و تتش کششی افزایش می‌یابد. اين امر منجر به 
اهر ای الم بسن نت بجر جهه فا خر ری نی 
شود. از سوی دیگر اين پوشش‌ها دارای غلظت بالایی 
از جای خالی یون اکسیژن هستند. که در دماهای بالا از 
یق تشکیل پوستهٌ اکسیدی محافظ پایدار (100- 
اکسیدی که در اثر حرارت رشد می کند) در سطح لایه 
وا نها اکسای با کستانیتر ی لاه #ابشط ۵ 
فص مرک ش شکب واسط سکس کید ابید 
امر منجر به جدایش و شکست ۲80ها خواهد شد 
[ ۱۶۰۱۵ ]. 
در پژوهش‌های انجام شده طی سالهای اخیر روش 
ها ملق برای هنم کام رزتیی او 
است [۱۶]. گزارش شده که پوشش‌های سد حرارتی بر 
پایه زیرکونات. به علت هدایت حرارتی پایین‌تر و 
پایداری فازی بالاتر و مقاومت بیشتر در برابر زینترینگ 
در دمای بالا گزینه مناسبی برای کاربردهای آینده در 
موتور هواپیماهاه توربین‌ها و دیگر قطعات دمای بالا 
باشند. (ضا) 02 حرضنا اخیرا به عتوان ماده:موزد استفاده 
در پوشش‌های سد حرارتی, پيشنهاد گردیده است 
[۱۵۰۱]. 
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با توسعه صنایع جدید. تقاضا برای پوشش‌های سد 
حرارتی با عملکرد عالی به یک ضرورت تبدیل شده 
است. در این زمینه» ارزیابی و پیش‌بینی طول عمر 
پوشش‌های سد حرارتی اهمیت بیشتری پافته است. 
عموماً طول عمر و دوام پوشش‌های 180 با تنش 
پسماند تعیین می‌شود. با توجه به آن که پوشش‌های سد 
حرارتی از چند لایه. شامل لایه فلزی و لایه سرامیکی؛ 
تشکیل شده‌اند تنش پسماند ناشی از عدم تطابق بین 
قراس انسیا راو فا اتکی شکنا. 
که پوشش تحت چرخه حرارتی مابین گاز داغ و هموای 
خنک‌کننده قرار می‌گیرد. اجتناب ناپذیر خواهد بود 
[۱۱۷۱۸] تنش پسماند منجر به جوانه‌زنی» رشد و 
گسترش ترک می‌شود [۸.۱۳,۱۹] و می‌تواند بر مدل 
شکشتت ,وشن سح خر ارشی: تسایر تک ازد راز اایتریه 
ارزیابی و پیش‌بینی الگوی توسعه تنش در پوشش سد 
حرارتی هنگام قرار گرفتن در معرض شوک حرارتی 
مهم و ضروری است [۱۰۱۷]. 

روش‌های شبیه‌سازی عددی به ویژه روش‌های 
المان محدود. در ارزیابی رفتار پوشش‌ها و لایه فصل 
مشترک تحت بارگذاری‌های حرارتی و مکانیکی از 
طریق محاسبه میدان‌های تنشی و کرنشی کاربرد فراوانی 
دارد و به منظور درک توزیع تتش درفصل مشترک فلز 
سرامیک مکانیزم گسترش تنش پسماند در لایه‌های 
مختلف و همچنین الگوی توسعه تسرک و فرآیند 
شکست 180 کاربرد فراوانی دارد [۱۱۷۰۱۹۰۲۵]. 

با توجه به مطالب ذکر شده. یافتن یک ترکیب 
جدید جهت جایگزینی با ۷۹2 به منظور بهبود خواص 
و دوام پوشش‌های سد حرارتی ضروری به نظر می- 
وسته لد ترش هایی در ناه را روا بهفتران 
کاندیدای جدید جهست ارتقای سطح کارایی این 
پوشش‌ها در کانون توجه پژوهشگران قرار گرفته اند. 
لیکن تاکنون پژوهش جامعی در ارتباط با پیش بینی 


رفتار حرارتی و توزیع تنش‌های پسماند در پوشش‌های 
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سد حرارتی بر پایه زیرکونات انجام نشده است. از 
اینرو در این تحقیق شبیه‌سازی عددی توزیع دماو 
تنش‌های پسماند در پوشش سد حرارتی سه لایه 
ترکیسی (192712007/8۷2/10۵1۷) تضت بسک 
چرخه حرارتی ترکیبی و واقعی. شامل گرمایش از دمای 
محیط تا دمای سرویس, نگهداری در دمای سرویس و 
سپس سرمایش تا دمای محیط مورد بررسی قرار گرفته 
است. 
مواد و روش تحقیق 

هندسه پوشش موردنظر متشکل از یک ساختار چندلایه 
است که از زیرلایه. لاه واسط لایه سرامیکی اولیه و 
نهایتاً لایه سرامیکی انوبه تضکیل شده امست, جنس 
لایههمای نام‌برده به ترتیب اینکونل ۷۳۸ 
ترکیب ۷0۸1۷ ترکیب سرامیکی زیرکونیای تقویت 
شده با ایتریا و نهایتاً ترکیب سرامیکی 122707 می- 
باشد که لایه واسط دارای ضخامت ۱۰۰ میکرومتر و 
سایر لایه‌ها دارای ضخامت ۲۰۰ میکرومتر هستند. 
نمایش گرافیکی شکل اصلی موتور توربوفن و مجموعه 
سیستم زیرلایه و پوشش, برش عرضی از مقطع دیواره 
و نمای سه بعدی مورد استفاده در شبیه‌سازی به ترتیب 
در شکل‌های (او ۲) آمده‌است. تمامی اجزای سیستم 
اعم از زیرلایه و پوشش‌هاء دوبعدی و شکل‌پذیر و با 
هندسه پوسته ای درنظر گرفته شده است. ضخامت لایه 
ها متناسب با نوع پوشش و نیز نوع لایه سطحی متفاوت 
است و از این رو در شکل (۲) بصورت شماتیک معین 
شده است. در این تحقیق. ضخامت هر لایه ۲۰۰ 


صفحه ۸-۸ 


زیر لایه 
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پوشش 


زبرلایه 


چ( 
شکل ۱ حوزه محاسبات انتخابی با توجه به هندسه واقعی جسم 
موردنظر 


برای ساده سازی محاسبات و جلوگیری از اعوجاج 
لخن از مه التاتها کل درک پم سای یهاش 
سرامیک /فلز و سرامیک/سرامیک تخت و بدون زبری 
فرض شده است. همچنین. فرض شده است که 
هیچگونه مقاومت حرارتی بین لایه‌ها وجود ندارد و از 


دانسیته تشثوری ترکیبات سرامیکی حین تحلیل استفاده 
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شده است (از تخلخل موجود در لایه‌های سرامیکی 
صرفنظر شده است). در تحلیل حاضر از المانهای مثلئی 
درجه دوم با شش گره در رئوس و وسط اضلاع استفاده 
شد. این نوع المان امکان بالاتر بردن صحت انتگرال- 
گیری برای محاسبات جابجاییهای مربوط به کرنش و 
تنش پسماند را فراهم می‌سازد. از طرف دیگر برای 
کاهش زمان حل از روش انتگرالگیری کاهشی از نوع 
کرنش صفحه ای بهره برده شده است. 

هنگام اعتصاص دادن حواص حرارتی و مکانیکی 
به لایه های مختلف. فرض همسانگردی خواص در هر 
دو بعد لحاظ شده است زیرا از حضور عیوب ناشی از 
فرایند پوشش‌دهی صرفنظر شده است. جدول (۱) 
حواص مورد استفاده در تحلیل کوپل حرارتی و 
مکانیکی مورد استفاده را بصورت تابعی گسسته از دما 
نشان می‌دهد. شایان ذکر است که نرم افزار اباکوس 
داده‌های بیشتر را با برونیابی استخراج و استفاده می‌کند. 
با توجه به بزرگ و قابل ملاحظه بودن حوزه محاسبات 
در قیاس با ابعاد عیوب بلوری حاضر در ساختار 
سرامیک‌ها؛ فرض همسانگردی لایه‌های پوشش قابل 


قبول است. 


شکل ۲ ترتیب لایه های پوشش سد حرارتی موردنظر در شبیه- 
سازی 
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جدول ۱ خواص فیزیکی و مکانیکی اجزای پوشش سد حرارتی 
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شکل ۲ محیط مش‌بندی‌شده مجموعه زیرلایه لایه واسط و لایه 


سرامیکی پوشش سد حرارتی 


سال بیست و هتم شماره دی ۱۳۹۲ 


شکل (۳) حوزه محاسبات را بصورت مش‌بندی 
شده نمایان می‌سازد. برای انجام تحلیل ساختاری و 
مکانیکی ضروری است. شرایط مرزی مختلفی بر روی 
نواحی کناری پوشش اعمال می‌گردد و در ایسن تحقیسق 
یزیر اس ابضرانوی خانپاه «سوایط مرزق مکتانک 
با اعمال قیود مناسب می‌توان از حرکت قطعه در جهات 
مقید شده جلوگیری نمود و نیز تفسیر صحبحی از 
تتشر‌هاین پمال درون ساده به دست آوود. از ایشرو: 
گوشه‌های افقی ساختار پوشش در جهت عمودی مقید 
شده و گوشه‌های عمودی از حرکت در جهت افقی 
بازداشته شده‌اند (شکل ۶). 

سیکل بارگذاری حرارتی بدین صورت است که 
ابتدا در مرحله گرم شدن مجموعه در مدت زمان ۳۰ 
ثانیه تا دمای مورد نظر سرویس گرم می‌شود و در 
مرحله سرویس به مدت یک ساعت در آن دما باقی 
می‌ماند. در آخرین مرحله با صرف زمان ۳۰۰ ثانیه 
سیستم حرارت خود را از دست خواهد داد. نمودار 
شماتیک بارگذاری حرارتی در شکل (۵) نشان داده 


شده است. 


۷ 


#9 


شکل ۶ اعمال شرایط مرزی مکانیکی به سیستم پوشش سد 


حرارتی 


نشریه مهندسی متالویژی و مواد 


زمان 


شکل ۵ بارگذاری حرارتی اعمالی روی سیستم پوشش سد حرارتی 


نتایج و بحث 
برای دستیابی به یک حل اقتصادی (با دقت مناسب و در 
کوتاهترین زمان ممکن) معمولاً می‌توان به روش‌های 
مختلف زمان لازم برای حل مساله را کاهش داد. یک 
کزشی ابو هو فافاشی ترش شاه الب هناب 
گونه‌ای است که حد پایداری افزایش یابد و تحلیل در 
تعداد گام‌های کمتری صورت پذیرد. این روش تکنیک 
مقیاس‌دهی به جرم (عصناههو دوعص نامیده می‌شود. در 
اکثر مساله ها از روش مقیاس دهی به جرم با مقیاس 
بالا برای کاهش زمان حل مساله بدون ایجاد خطا در 
جواب استفاده می‌شود. مقایسه نمودارهای انرژی داخلی 
و جنبشی معیاری برای ارزیابی عدم ایجاد خطا در 
مساله با استفاده از این روش است. میزان گام زمانی 
پایدار در نرم افزار به وسیله تخمین آن بر روی تک تک 
المان‌ها و سپس انتخاب کوچک ترین گام زمانی تعیین 
می‌گردد. به همین دلیل در صورتی که در مدل از المان- 
های بسیار کوچک استفاده شده باشد یابه هر دلیل 
قسمتی از مدل به دلیل تغییر شکل های پیچیده طول 
مشخصه کوچکی را به حل تحمیل کند. مقدار گام زمانی 
بسیار کوچک می گردد و در نتیجه تعداد گام زمانی مورد 
نیاز برای کل زمان شبیه سازی و به دنبال آن زمان صرف 
شده برای شبیه سازی بسیار طولانی می شود. گاهی 
آوقات این المان های کوک تعدادشان بسیان آندک بودة 
و در قسمت کوچکی از کل مدل متمرکز هستند به 
صورتی که تاثیر چندانی در دقت حل نداشته ولی به 
شدت باعث کم شدن حد پایداری می‌گردند. در این 


نشريه مهندسی متالورژی و مواد 


شبیه سازی و بهینه سازی حل عددی... 


شرایط لازم است با استفاده از روش‌هایی این گام زمانی 

پایدار را افزایش داد. به گونه ای که علاوه بر زمان حل؛ 

دقت مساله نیز تحت تاثیر قرار نگیرد. با توجه به 

معادلات (۱) و (۲), مشخص می‌گردد که افزایش 
۳3 


162 (۱) 


۳( 5 2 0 
با توجه به فرمولهای بالاء با افزایش میزان چگالی 
ماده به میزان مورد نظر» سرعت حرکت صوت درون 
جسم به میزان جذر آن مقدار کاهش می‌یابد و به همان 
میزان نمو حلی افزايش می‌یابد. این افزایش نموهای 
حلیء کل زمان حل را به بازه های زمانی کوچکتری 
تقسیم می‌کند. حلگر ضمنی نرم‌افزار با استفاده از همین 
اصل؛ روشی را برای کاهش زمان شبیه‌سازی ارائه می- 
دهد که به مقیاس‌دهی به جرم معروف است. در این 
روش به طور مجازی. جرم کل مدل یا قسمتی از آن 
افزایش يافته و به این ترتیب گام زمانی پایدار افزایش می 
یابد. در صورتی که مقیاس دهی به جرم به درستی انجام 
شود می‌تواند علاوه بر این که زمان حل را کاهش دهد 
دقت حل را نیز در حد قابل قبولی حفظ کند. در صورتی 
که شبیه‌سازی شبه استاتیک باشد. کار انجام شده توسط 
نیروهای خارجی تقریباً برابر با انرژی داخلی سیستم 
خواهد بود. در اپن حالت انرژی ویسکوز تلف شده 
معمولاً کوچک است. مگراینکه از مواد ویسکو الاستیک 
استفاده شده باشد. همانطور که از ماهیت فرآیند های 
شبه استاتیک انتظار می رود به دلیل سرعت بسیار کم 
ماده درون مدل. نیروهای اینرسی بسیار کوچک و قابل 
صرف نظر کردن هستند . به همین دلیل انرژی جنبشی 
ماده قابل تغییر فرم نیز بسیار کوچک خواهد بود .لذا 
طبق یک قانون عمومی. انرژی جنبشی ماده تغییر شکل 
پذیر نباید از بخش کوچکی (معمولاً ۵ تا۱۰ درصد) از از 
انرژی داخلی آن در حین فرآیند تجاوز کند. 
حین حل عددی از قانون کلی ذکر شده برای 
بررسی صحت شبیه‌سازی شبه استاتیک استفاده می‌شود. 
در ابتدای تحلیل از یک مقیاس بندی جرم نسبتاً بزرگ 


استفاده شده و شبه‌سازی انجام می‌شود. پس از دستیابی 


سال بیست و هشتم شماره یک؛ ۱۳۹۵ 
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به نتایج بایستی نسبت انرژی جنبشی به انرژی داخلی 
ماده رسم و بررسی شود . در صورتی که نسبت این دو 
نرژی از مقدار گفته شده تجاوز نکرده باشد. می‌توان 
طمینان حاصل نمود که نتایج مربوط به یک فرآیندشبه 
ستاتیک است. در غیر این صورت مقدار مقیاس دهی به 
جرم کاهش داده می شود و دوباره مساله حل می گردد و 
ین فرآیند تا جایی ادامه می یابد که نسبت مذکور در 
بیشتر زمان فرایند به حد مورد نظر برسد و در عین حال 


ختلاف پاسخ‌ها بین دو شبیه‌سازی متوالی ناچیز گردد. 


فزایش انرژی جنبشی به مقداری بیش از مقدار ذکر شده 
بیانگر این است که مقداری از انرژی صرف شده برای 
غلبه به اینرسی ماده و در نتیجه شتاب دادن به جسم 
صرف شده است. شکل (1) نمودارهای تغییرات انرژی 
جنبشی و داخلی را در محدوده معینی از فرآیند که 
تغییرات دما به حد پایا می‌رسد. نشان می‌دهد. همانطور 
که در این شکل نشان داده شده است. انرژی جنبشی 
کسر بسیار اندکی از کل انرژی درون سیستم را صرف 
می‌کند و حالت شبه استاتیک سیستم قابل دسترسی 
است. روند تغییرات در ۰ شثانیه حرارتدهی رسم شده 
است و این روند در بقیه فرایند ادامه دارد که در این 
شگل برای نشان دادن بهتر مقایسه. فقط همان مرحله 


اول نشان داده شده است. 


و و۱۱ 


و ۱۰6۱ 


تلا 


اعوجاج مش‌ها در مدل کوپله المان محدود تا 
جایی ادامه می‌یابدکه مش از به دست اوردن نتایج دفیق 
عاجز است و در نتیجه تحلیل به دلایل زیادی خاتمه می- 
پابد .در چنتین مسائلی ضروری است که از ابزار سش 
بندی هوشمند و تغییریافته‌ای برای کاهش موثر 
تغییرشکل و اعوجاج مش‌ها و مشکلات ناشی از آن 
ستفاده گردد. مش‌بندی تطبیق‌پذیر ایزاری است که بساقی 
ماندن مش با کیفیت را حتی با وجود تغییر فرم‌های شدید 
در طول تحلیل امکان پذیر می‌سازد. نرم‌اف زار المان 
محدود قابلیت‌های بسیار زیاد و قوی مش‌بندی تطبیق 
پذیر را برای تحلیل مسائل به شدت غیرخطی, در دامنه 
مسائل شبه استاتیکی تا مسائل دینامیک فراهم کرده 
است. در این تحقیق از تنظیم دو رابطه فرکانس و میزان 
تغییر گره‌ها در هر نمو حلی برای حداقل نمودن میزان 
اعوجاج مش‌بندی و همچنین حداقل نمودن زمان حل 
مدل استفاده شد. 

قر کاس بیانگر این سباله هس که بعد از ین 
نمو حلی گره‌ها به مکان اصلی بازگردادنیده شوند.این 
امر با امکان دادن به مش برای حرکت نسبتاً مستقل از 
حرکت ماده انجام می پذیرد که باعث می‌گردد مش با 
هندسه جدید قطعه کار تطبیق یابد با این وجود مش- 
بندی تطبیق‌پذیر تعداد المان‌ها و اتصال بین آن‌ها را تغییر 


نمی‌دهد. 
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شک 1 مقایسه تغییرات انرژی جنبشی و داخلی حین فرایند 
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(ب) 


شکل ۸ کانتورهای الف) توزیع تنش پسماند؛ و ب) اعوجاج لایه‌ای در سیستم پوشش سد حرارتی 


بعد از چندین بار تکرار حل مساله. کمترین میزان 
تغییرات هندسه مش و بیشترین میزان کاهش زمان حل 
نآزا فرکانین ۹9۱ (تخاک ففعانت نمی گرمهر 
هر فرکانس) بدست آمد. به بیان دیگرء با رعایت 
معیارهای همگرایی در محدوده‌ای از حل از تحلیل 
لاگرانژی استفاده شده و هندسه ماده به مش‌بندی مقید 
است و در ادامه حل. دیگر ماده و مش‌بندی مستقل از 
یکدیگر خواهند شد. از جایی به بعد با ادامه حل؛ 
فضای اویلری حین پیشرفت مساله حکفرم ا است و 


نشریه مهندسی متالورژی و مواد 


احتمال صفرشدن جرم المان را از بين می‌برد. 

شکل (۷) توزیع حرارت و شار حرارتی جاری 
درون لایه‌ها و شکل (۸) توزیع تنش پسماند و میزان 
اعوجاج الاستیک درون مجموعه سیستم پوشش سد 
بفزازتی زا نانمشد 

همانطور که در شکل‌های (۷-الف) و (۷-ب) دیده 
مي‌ شسود پیشترین مان مر که حرارتی در در لاه 
پوشش‌های سد حرارتی فوقانی دیده می‌شود که با 
فلسفه طراحی و زمینه کاربرد این دسته از سیستم‌های 
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محافظ در توافق است. پوششهای سرامیکی دو لایه با 
دمای ذوب بسیار بالا بعنوان سد حرارتی عمل می‌کنند 
و بیشترین میزان انرژی حرارتی را درون خود حفظ 
می‌کنند تا به لایه‌های زیرین که از جنس فلزی هستند و 
دماهای نرم شدن, خمیری شدن و پارگی داغ کمتری 
دارند. توا کمتری وارد گردد [۲1]. بعلاوه همانطور 
که در شکل (۷- ب) دیده می‌شود زیرلایه از جنس 
سوپرآلیاژ نیز متحمل کمترین میزان بار حرارتی می‌گردد 
و حداکثر عمر را برای توربوفن به همراه خواهد داشت. 

در شکل (۸) کارکرد مثبت دیگر این گونه طراحی 
پوشش‌های ترکیبی سد حرارتی دیده می‌شود. حداکثر 
مقادیر تنش پسماند در درون محدوده دو پوشش 
سرامیکی بالایی قرار گرفته است و تقریبا تنش‌های 
مخرب کاهنده عمر که ناشی از تغییرات دمایی در 
میقم هسفق ان لانه واسط بو تمی کش و این اسر 
عمر زیرلایه را افزايش می‌دهد. با مقایسه تحلیل‌های 
مشابه که روی پوشش‌های سد حرارتی ساده تک لایه 
انجام شده است [۲۱] می‌توان به اين نتیجه رسید که 


میزان تنش ذخیره شده در زیرلایه با سیستم‌های تک 


1 
030 940 050 


لایه بیشتر از سیتم‌های نوین ارائه شده در این تحقیق 
است. باید به مزایای فوق این مزیت را نیز اضافه نمود 
که تمرکز شارحرارتی در لایه‌های بالایی باعث حصول 
زودهنگام تعادل حرارتی در نیمه بالایی سیستم پوشش 
شده و کاهش حجم حرارت انتقالی به ترکیبات فلزی 
زیرین را به همراه خواهد داشت. 

شکل )٩(‏ میزان اعوجاج در راستای شعاعی را 
برای مجموعه سیستم زیرلایه و سه پوشش بالایی 
متشکل از یک لایه واسط و دو لایه سرامیکی سد 
حرارتی نشان می‌دهد. همانگونه که در این شکل دیده 
مي‌شود؛ تمامی افت دمایی و انرژی حرارتی بصورت 
قابل ملاحظه‌ای در دو لایه سرامیکی بالایی متمرکز شده 
است و حتی در لایه اتصال هم چندان بار حرارتی و 
تخریب ناشی از آن را مشاهده نمی‌کنيم. مقایسه این 
نتایج با نتایج حاصل از شبیه‌سازی پوششهای سد 
حرارتی معمولی. افت دمایی بیشتر و دمای زیرلایه 
پایین‌تر را برای این پوششهای نوین نمایان می‌سازد 
[ ۱۶۰۱۵ ]. 


].2-6[ 


اعوجاج الاستیک (میکرومتر ) 
۳ 
3 
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فاصله از سطح پوشش یه سمت زیرلایه (میکرومتر) 


شکل ٩‏ الگوی توزیع اعوجاج در طول سیستم پوشش و لایه‌های پوشش سد حرارتی 
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افزوده شدن ضخامت به زبرلایه و پوشش از سطح زیرلابه (میلیمتر) 


شکل ۰ روند افزایش دمای سطح سیستم از زیرلایه به سمت پوشش‌ها 


شکل (۱۰) روند افزایشی تغییرات دمایی را از 
همانطور که در این شکل به وضوح مشاهده می‌شود در 
زیرلایه از جنس سوپرآلیاژ شاهد حداقل دماهای کاری 
در سیستم هستیم که این خود به افزايش عمر توربوفن 
کمک می‌کند. در ادامه با واردشدن به ناحیه لایه واسط 
و لایه‌های دوتایی سرامیکی بعدی, دو اتفاق مثبت در 
راستای ازدیاد بیشتر عمر زیرلایه فلزی می‌افتد. نخست 
اینکه. در لایه‌های فوقانی پوشش شاهد دماهای بالاتر 
تا زیرلایه قابل ملاحظه می‌باشد. این امر به نوبه خحود 
کارکرد مثبت حضور سیستم‌های سد حرارتی را در 
کاهش دمای زیرلایه نشان می‌دهد. مزیت دوم تمرکز 
بیشترین میزان بار حرارتی در محدوده ضخامت دو لابه 
که این امر در مقایسه با پوشش سد حرارتی تک لایه 
سرامیکی نشان دهنده قابلیت عایق‌سازی قابل توجه 
پوشش دولایه ترکیبی پيشنهاد شده در پژوهش حاضر 
لایه نشانی» علاوه بر نشان دادن مقاومت به زینتر بالاتر 


نشريه مهندسی متالورژی و مواد 


و از دست نرفتن خواص عایق حرارتی در دمای بالاء به 
دلیل دولایه بودن و چسبندگی مناسب. می‌تواند مقاومت 


بیشتری در برابر بار حرارتی از خود نشان دهد [۱ ۲]. 


نتیجه گیری 
در این تحقیق شبیه‌سازی عددی توزیع دما و تنش‌های 
پسماند در پوشسش سد حرارتسی سسه لاه 
 ۷(‏ تحت یک جرخه 
تا دهای سرویس» نگهدذ اوق دز همای سرویین و نسیسن 
کانتورهای دمایی و شار حرارتی نشان داد که پوششهای 
گرمای کمتری را به زیرلایه انتقال می‌دهند و از انتقال 
انرژی حرارتی نیز جلوگیری می‌کنند. توزیع تنش 
پسماند در طول منطقه پوشش-زیرلایه بصورت عمده 
در ناحیه پوششهای سرامیکی متمرکز شده است و 
بیشترین میزان اعوجاج نیز در همین ناحیه قرار دارد. 
این مطالب گویای این است که پوشش‌های دو لایه سد 


حرارتی دارای بیشترین میزان کاربرد در افزایش عمر 


سال بیست و هشتم شماره یک؛ ۱۳۹۵ 


نسیم نایب پاشایی - حسین وفایی نژاد- سید حسین سیدین- محمدرضا ابوطالبی - سید محمد مهدی هادوی- حسین سرپولکی ۱۱ 


قطعات توربین ها هستند. بعلاوه از تکنیک افزایش . تطبیق‌پذیر انتقال مش از حالت لاگرانژی به اویلری از 
مجازی جرم و دانسیته المان برای افزایش نمو حل و اعوجاج بیش از حد و صفرشدن جرم المان جلوگیری 
تیان آزا کاملی کل وان سخایشا تیا مرفتیت بهسره ...کرد و ریش گزآبسامساله و تفا کریتل با 
برده شد. علاوه بر آن. استفاده از سیستم هوشمند و اطمینان بیشتری حاصل شد. 
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